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Spins  confined  in  quantum  dots  are  a  leading  candidate  for  solid  state  quantum  bits  that  can 
be  coherently  controlled  by  optical  pulses.  There  are,  however,  many  challenges  to  developing 
a  scalable  multi-bit  information  processing  device  based  on  spins  in  quantum  dots,  including  the 
natural  inhomogeneous  distribution  of  quantum  dot  energy  levels,  the  difficulty  of  creating  all-optical 
spin  manipulation  protocols  compatible  with  nondestructive  readout,  and  the  substantial  electron- 
nuclear  hyperfine  interaction-induced  decoherence.  Here,  we  present  a  scalable  qubit  design  and 
device  architecture  based  on  the  spin  states  of  single  holes  confined  in  a  quantum  dot  molecule.  The 
quantum  dot  molecule  qubit  enables  a  new  strategy  for  optical  coherent  control  with  dramatically 
enhanced  wavelength  tunability.  The  use  of  hole  spins  allows  the  suppression  of  decoherence  via 
hyperfine  interactions  and  enables  coherent  spin  rotations  using  Raman  transitions  mediated  by  a 
hole-spin-mixed  optically  excited  state.  Because  the  spin  mixing  is  present  only  in  the  optically 
excited  state,  dephasing  and  decoherence  are  strongly  suppressed  in  the  ground  states  that  define 
the  qubits  and  nondestructive  readout  is  possible.  We  present  the  qubit  and  device  designs  and 
analyze  the  wavelength  tunability  and  fidelity  of  gate  operations  that  can  be  implemented  using 
this  strategy.  We  then  present  experimental  and  theoretical  progress  toward  implementing  this 
design. 

PACS  numbers:  78.20.Ls,  78.47.-p,  78.55.Cr,  78.67.Hc 


Single  confined  spins  have  long  been  considered  as 
possible  bit  states  for  novel  optoelectronic  logic  devices, 
including  quantum  computers—  Spins  confined  in  III-V 
semiconductor  self-assembled  quantum  dots  (QDs)  have 
received  a  great  deal  of  attention  because  they  interact 
strongly  with  light  and  provide  the  opportunity  for  ul¬ 
trafast  all-optical  implementation  of  logic  operations— 
There  has  been  dramatic  progress  in  the  initialization,  co¬ 
herent  manipulation  and  readout  of  single  spins  in  GaAs 
and  InGaAs  QDs*£A2  but  many  challenges  to  the  cre¬ 
ation  of  a  scalable  quantum  logic  device  based  on  optical 
control  of  single  spins  remain.  One  serious  obstacle  is  the 
natural  inhomogeneous  distribution  of  energy  levels  in  a 
quantum  dot  ensemble.  This  distribution  necessitates  in¬ 
dividually  tuned  lasers  to  control  each  bit  and  restricts 
the  ability  to  mediate  entanglement  via  photonic  cavity 
modes— -  Another  obstacle  is  that  existing  approaches  to 
all-optical  coherent  control  of  spins  require  a  transverse 
magnetic  field  in  order  to  mix  spin  states  and  permit 
Raman  transitions—  The  transverse  magnetic  field  pre¬ 
vents  nondestructive  readout—  A  third  obstacle  is  the 
suppression  of  decoherence.  Among  the  many  channels 
for  decoherence,  the  hyperfine  interaction  with  nuclei  is 
the  most  detrimental—^-17 

We  present  a  quantum  bit  (qubit)  design  and  device 
architecture  that  overcomes  many  of  these  obstacles.  In 
our  approach,  the  two  spin  states  of  a  single  hole  en¬ 
code  the  qubit.  The  qubit  hole  is  localized  in  the  top 
QD  of  a  coupled  pair  of  vertically-stacked  InAs  QDs. 
These  stacked  pairs  of  QDs  are  known  as  Quantum  Dot 
Molecules  (QDMs)  because  coherent  tunneling  leads  to 


the  formation  of  states  with  delocalized  wavefunctions 
that  have  molecular  symmetries , 18-21  We  propose  initial¬ 
izing  and  rotating  the  spin  using  optically  excited  delo¬ 
calized  molecular  states  with  mixed  spin  character.  A  dif¬ 
ferent  optically  excited  state  that  is  immune  to  spin  mix¬ 
ing  provides  a  recycling  transition  for  readout.  We  show 
that  this  approach  enables  all-optical  coherent  control 
of  single  hole  spins,  controlled  interactions  between  two 
non-identical  qubits  mediated  by  photonic  cavity  modes, 
and  nondestructive  readout  of  the  spin  projection  along 
the  growth  axis. 

The  approach  presented  here  builds  on  experimen¬ 
tal  progress  in  the  coherent  control  of  single  hole 
spins16i17i22~-^  and  electron  spins  in  QDMs—  The  new  de¬ 
sign  combines  advantages  of  these  approaches  and  makes 
several  key  improvements  that  significantly  enhance  scal¬ 
ability.  First,  the  hole  spin  is  less  sensitive  to  dephasing 
by  hyperfine  interactions  with  nuclei,  which  provide  a 
dominant  decoherence  mechanism  for  electrons:  16i28  Sec¬ 
ond,  the  spin  mixing  for  hole  tunneling  in  QDMs  can 
be  engineered  to  be  much  larger  than  the  mixing  for 
electron  tunneling,  greatly  enhancing  the  fidelity  of  the 
optically-driven  qubit  rotation.  Third,  the  spin  mixing 
in  QDMs  provides  a  pathway  to  implement  all-optical 
coherent  control  using  only  magnetic  fields  applied  along 
the  growth  direction,  which  is  parallel  to  the  optical  axis. 
As  a  result,  the  approach  is  compatible  with  the  use 
of  recycling  transitions  for  nondestructive  readout.  Fi¬ 
nally,  our  approach  utilizes  optical  transitions  between 
electrons  and  holes  located  in  separate  quantum  dots. 
These  transitions  have  an  extremely  large  Stark  shift  as 
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a  function  of  applied  electric  fields,  and  thus  provide  an 
enhanced  capacity  to  tune  the  transitions  into  resonance 
with  external  laser  sources  or  optical  cavity  modes. 

The  use  of  ‘indirect’  optical  transitions,  involving  elec¬ 
trons  and  holes  in  separate  QDs,  underlies  an  important 
conceptual  aspect  of  the  device  architecture  proposed 
here.  Existing  device  designs  often  assume  an  ensemble 
of  identical  QDs,  which  is  impossible  to  achieve  experi¬ 
mentally.  Demonstrations  of  spin  initialization,  control, 
or  readout  using  QDs  typically  identify  the  specific  opti¬ 
cal  transitions  of  a  particular  QD  and  then  tune  optical 
cavities  and  external  laser  sources  into  resonance  with 
that  particular  QD.  This  approach  cannot  be  scaled  be¬ 
yond  a  small  number  of  qubits.  In  our  approach  we  begin 
by  accepting  the  inhomogeneous  distribution  of  QD  en¬ 
ergy  levels.  We  design  a  device  architecture  in  which 
the  optical  transitions  of  a  subset  of  individual  QDMs 
within  this  ensemble  can  be  tuned  into  resonance  with  a 
fixed-frequency  optical  cavity  or  external  laser.  The  key 
to  tuning  many  individual  QDMs  into  resonance  with  a 
fixed-frequency  optical  cavity  or  laser  is  the  use  of  in¬ 
direct  transitions  whose  energy  shifts  strongly  with  the 
electric  field  applied  locally  to  individual  QDMs. 

We  present  detailed  modeling  of  the  QDM  states  that 
could  be  used  for  this  scheme,  compute  the  fidelity  of  gate 
operations  implemented  with  this  wavelength-tunable 
approach  and  provide  experimental  evidence  of  the  large 
spin  mixing  necessary  for  such  a  scheme.  In  Sect.  Q]  we 
present  our  qubit  design  and  describe  how  this  design 
suppresses  hole  spin  decoherence  or  dephasing  in  the  log¬ 
ical  basis  states.  In  Sect.  EH  we  develop  our  strategies  for 
full  single-qubit  all-optical  control.  We  first  analyze  the 
use  of  indirect  optical  transitions  (subsection  III  All  and 
then  describe  how  these  indirect  transitions  enable  single 
spin  initialization  and  readout  with  enhanced  wavelength 
tunability  (subsectionlll  Bl).  In  subsectionlll  Clwe  develop 
the  coherent  control  of  the  hole  spin  using  Raman  transi¬ 
tions  mediated  by  optically  excited  states  with  hole  spin 
mixing.  In  Sect.  IIIII  we  address  scalability  in  the  con¬ 
text  of  wavelength  tunability  and  fidelity  of  spin  control 
protocols  and  discuss  how  the  conditional  interactions 
can  be  mediated  by  photonic  cavity  modes.  In  Sects.  IIVI 
and  [V]  we  present  progress  toward  the  implementation  of 
this  device  design,  including  experimental  measurements 
of  hole-spin  mixing  in  the  optically  excited  positive  trion 
state  and  theoretical  calculations  of  QDM  structures  that 
enhance  the  spin  mixing. 


I.  QUBIT  DESIGN 

In  Fig.  [ljr  we  schematically  depict  the  band  structure 
and  lowest  single  particle  energy  levels  for  electrons  and 
holes  in  the  QDM  designed  for  this  application.  The 
QDM  grown  by  Molecular  Beam  Epitaxy  is  embedded 
within  a  p-i-n  diode  structure  that  allows  an  applied  elec¬ 
tric  field  to  tune  the  energy  levels  of  the  QDs  relative  to 
one  another  and  the  Fermi  level.  The  bottom  QD  is 


truncated  to  a  smaller  height  than  the  top  QD  so  that 
the  hole  levels  can  be  tuned  into  resonance  while  elec¬ 
tron  levels  retain  a  significant  energy  offset  that  causes 
electrons  to  always  relax  to  the  top  QD^  The  tunnel 
barrier  formed  by  the  intrinsic  GaAs  region  between  the 
p-doped  substrate  and  the  bottom  QD  is  chosen  so  that 
the  QDM  remains  deterministically  charged  with  a  sin¬ 
gle  hole  in  the  electric  field  range  of  operation.  Tunneling 
of  electrons  into  the  QDM  is  prevented  by  the  inclusion 
of  an  ALcGa^-^) As  ( x  <  0.4)  blocking  layer.  The  het¬ 
erostructure  design  is  schematically  depicted  in  Fig.  []Jl 
and  numerical  values  for  the  design  are  discussed  in  Ap¬ 
pendix  [A] 
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FIG.  1:  (Color  Online)  a)  Schematic  band  diagram  of  pro¬ 
posed  qubit  QDM.  Direct  (D)  and  Indirect  (I)  optical  transi¬ 
tions  are  schematically  indicated,  b)  Cross-sectional  STM  of 
a  QDM  demonstrating  the  stacking  of  QDs  and  the  possibil¬ 
ity  of  lateral  offsets,  c)  Photonic  cavity  architecture  incorpo¬ 
rating  multiple  QDMs.  d)  Schematic  cross  section  of  device 
heterostructure . 


Our  qubit  is  defined  by  the  spin  projections  along  the 
growth  axis  of  a  single  hole  (h+)  confined  in  the  QDM. 
The  optically  excited  state  has  one  additional  electron- 
hole  pair  and  is  called  a  positive  trion  (X+).  We  use 
ihVhr)  to  desci'ibe  the  spatial  location  and  spin  orien¬ 
tation  of  each  charge:  es  ifir)  are  the  electron  spin 
projections  in  the  bottom  (top)  QD;  j\  i  correspond  to 
Se  =  ±1/2.  Holes  in  QDs  contain  both  light-  and  heavy- 
hole  contributions  and  hole  spins  are  properly  described 
as  Luttinger  spinors^  In  single  InAs  QDs,  the  heavy- 
hole-only  approximation  is  largely  valid,  though  the  con¬ 
tribution  of  light  holes  does  impact  spin  dynamics  In 
QDMs  the  contribution  of  light-hole  states  becomes  more 
important  and  leads  to  unique  and  tunable  properties 
for  hole  spins  i^2r—  The  hole  spin  mixing  that  emerges 
in  QDMs  is  a  key  element  of  the  approach  we  present 
here^Si33,35  Although  the  hole  spinors  contain  contribu¬ 
tions  from  all  light-  and  heavy-hole  spin  projections,  the 
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spinors  are  dominated  by  a  single  heavy-hole  spin  pro¬ 
jection.  For  clarity  we  label  the  hole  state  as  ft,  (l,  cor¬ 
responding  to  the  dominant  heavy  hole  spin  projection 
(Jz  =  ±3/2)  in  each  QD.  The  (////)  notation  describes 
the  states  of  the  QDM  far  away  from  resonances.  When 
the  electric  field  tunes  energy  levels  into  resonance  the 
resulting  molecular  states  can  be  described  as  symmetric 
and  antisymmetric  combinations  of  these  basis  states. 

In  Fig.  [5Ji  and  b  we  plot  the  calculated  energy  levels  of 
the  h+  and  X+  states  as  a  function  of  the  applied  electric 
field  in  the  presence  of  a  static  1  T  magnetic  field  applied 
along  the  growth  axis  and  optical  axis  (Faraday  geome¬ 
try).  The  energy  level  calculation  uses  matrix  Hamilto¬ 
nian  methods  that  have  been  shown  to  accurately  model 
the  states  of  QDMs^21i32'35  —  The  matrix  Hamiltoni¬ 
ans  and  numerical  values  used  in  the  calculations  pre¬ 
sented  here  are  described  in  Appendix  [B]  We  focus  first 
on  the  energy  levels  for  the  hole  (h+).  There  are  four 

hole  states:  We  reference  en¬ 

ergies  to  the  energy  of  a  hole  in  the  top  QD,  so  the  two 
states  with  a  hole  in  the  top  QD  have  an  energy  that  does 
not  depend  on  electric  field  (horizontal  lines  in  Fig.  [2}r). 
States  with  a  hole  in  the  bottom  QD  have  a  linear  de¬ 
pendence  of  energy  on  the  applied  electric  field  (diagonal 
lines  in  Fig.  [2]d) .  When  the  hole  levels  are  in  resonance, 
coherent  tunneling  leads  to  the  formation  of  molecular 
orbitals  and  the  appearance  of  anticrossings,  as  shown  in 
Fig.  dr. 

The  formation  of  molecular  orbitals  results  in  pertur¬ 
bations  to  the  hole  spin  g  factor  and  creates  the  elec¬ 
tric  field  dependent  Zeeman  splitting  that  can  be  seen 
in  Fig.  O 3*22  In  the  absence  of  hole  spin  mixing,  the 
coherent  tunneling  that  results  in  molecular  orbitals  is 
spin-conserving  and  the  four  hole  states  retain  their  spin 
values  throughout  the  range  of  electric  field.  The  pres¬ 
ence  of  a  lateral  offset  between  QDs  (Fig.  [TJd)  ,  however, 
breaks  the  QDM  symmetry  and  creates  an  effective  spin- 
flip-tunncling  mechanism  that  mixes  hole  states  with  op¬ 
posite  spin  projections  located  in  separate  QDs<22  The 
magnitude  of  hole  spin  mixing  included  in  these  calcula¬ 
tions  (see  Appendix  [Bj  does  not  result  in  the  appearance 
of  anticrossings  in  Fig.  because  the  magnetic  field  is 
too  small  to  bring  the  spin  mixed  states  sufficiently  close 
in  energy.  The  presence  of  the  spin  mixing,  however, 
appears  in  calculations  of  the  hole  eigenstate  composi¬ 
tion.  In  Fig.  [2jc  we  plot  the  relative  contribution  of  each 
of  the  four  hole  basis  states  to  the  molecular  hole  state 
dominated  by  ''j .  Within  certain  ranges  of  electric 
field,  the  hole  spin  mixing  leads  to  eigenstates  containing 
significant  contributions  from  multiple  spatial  and  spin 
configurations. 

We  design  our  qubit  and  device  architecture  to  take 
advantage  of  hole  spin  mixing  in  the  optically  excited 
state  for  initialization  and  control  while  suppressing  the 
potential  negative  consequences  of  hole  spin  mixing  in 
the  logical  basis  states  where  information  is  stored.  The 
logical  ground  states  are  the  spin  projections  of  a  sin- 
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FIG.  2:  (Color  Online)  a,b)  Calculated  energy  levels  for  the 
X+  (a)  and  h+  (b)  states  as  a  function  of  applied  electric 
field  in  the  presence  of  a  1  T  magnetic  field  parallel  to  the 
optical  axis.  Right  column  indicates  the  dominant  atomic- like 
basis  states  for  each  group  of  molecular  states,  c)  Hole  spin 
contributions  to  one  hole  eigenstate  as  a  function  of  applied 
electric  field. 


gle  hole  confined  in  a  single  QD.  The  use  of  holes  sup¬ 
presses  hyperfine  interactions  with  nuclei  and  leads  to 
dephasing  times  estimated  to  be  at  least  100  ns  for  a 
hole  in  a  single  QD*22  For  electric  fields  to  the  left  of  Fp 
in  Fig.  [2J  the  logical  basis  states  approach  the  electric 
field  of  coherent  tunneling.  The  formation  of  molecular 
states  may  have  negative  consequences  for  the  storage  of 
quantum  information,  including  increased  spin-orbit  in¬ 
teraction,  electric-field  induced  changes  in  g  factor  and 
hole  spin  mixing42  To  minimize  the  potential  negative 
impact  of  these  effects  on  the  storage  of  quantum  infor¬ 
mation,  we  design  the  device  architecture  to  operate  only 
in  the  range  of  electric  fields  to  the  right  of  Fp  in  Fig.  [2] 
where  the  logical  ground  states  are  localized  atomic-like 
states  with  minimal  perturbations  of  their  spin-orbit  in¬ 
teraction,  g  factor  or  hole  spin  mixing^  In  Fig.  [5J:  we 
show  that  the  logical  basis  states  remain  pure  for  electric 
fields  to  the  right  of  Fp:  the  molecular  state  is  dominated 

by  with  contributions  from  (^’o)  and  below 

1%  and  the  contribution  of  below  0.1%. 

The  addition  of  an  electron-hole  pair  to  the  QDM 
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changes  the  Coulomb  interactions  and  thus  the  electric 
field  at  which  coherent  tunneling  leads  to  the  formation  of 
molecular  states.  Our  proposed  all-optical  control  strat¬ 
egy  operates  in  the  range  of  electric  fields  where  the  logi¬ 
cal  basis  states  remain  sufficiently  pure  but  the  optically 
excited  state  can  take  advantage  of  hole  spin  mixing.  As 
we  describe  below,  hole  spin  mixing  in  the  optically  ex¬ 
cited  state  allows  us  to  perform  all  of  our  spin  opera¬ 
tions  without  a  transverse  magnetic  field.  Consequently, 
spin  projections  along  the  optical  axis  are  eigenstates 
of  the  Hamiltonian  and  our  design  is  compatible  with 
nondestructive  readout.  In  the  next  sections  we  develop 
the  strategy  for  spin  initialization,  coherent  control,  and 
readout. 


II.  FULL  SINGLE  QUBIT  OPTICAL  CONTROL 

We  propose  to  manipulate  the  qubit  states  with  op¬ 
tical  pulses  that  couple  to  trion  (X+)  states.  The  X+ 
state  contains  one  electron  located  in  the  top  QD  and 
two  holes  that  can  be  in  either  the  top  or  bottom  QD. 
Coulomb  interactions  with  the  additional  electron  and 
hole  cause  the  anticrossings  for  the  X+  states  to  happen 
in  a  different  range  of  electric  fields  than  the  anticross¬ 
ings  of  the  h+  states,  as  shown  in  Fig.  [2]  There  are  12 
basis  states  for  the  X+  with  unique  spatial  and  spin  dis¬ 
tributions.  The  two  states  that  have  two  holes  in  the 
bottom  QD  are  outside  the  energy  range  we  must  con¬ 
sider,  so  only  10  states  appear  in  Fig.  [2ji,  as  described 
in  Appendix  iBl  If  both  holes  are  located  in  the  top  QD, 
the  Pauli  exclusion  principle  requires  them  to  be  in  a 
spin  singlet,  e.g.  .  If  the  two  holes  are  in  separate 

QDs,  however,  both  singlet  and  triplet  configurations  are 
possible.19!21  We  denote  the  singlet  state  as  and 

the  three  triplet  states  as  and  An 

analogous  set  of  states  exist  for  the  electron  spin  down 
case  (j,). 

Coherent  coupling  of  any  two  states  leads  to  anticross¬ 
ings  of  the  energy  levels.  If  only  spin-conserving  tun¬ 
neling  was  possible,  only  the  singlet  states  would  tunnel 

couple  (e.g.  (0’^)  -H-  )  and  triplet  states  would 

pass  through  the  resonance  without  coupling.  In  Sect.HVl 
we  show  experimental  evidence  for  the  existence  of  hole 
spin  mixing  that  couples  singlet  and  triplet  states.  This 
hole  spin  mixing  is  included  in  the  calculated  energy  lev¬ 
els  presented  in  Fig.  [2]  Asa  result  of  this  mixing,  many  of 
the  trion  states  that  appear  in  Fig.  [2}i  are  molecular-like 
admixtures  of  several  basis  states.  We  take  advantage  of 
this  mixing  to  enable  new  optical  control  strategies. 

Figure  0]  presents  a  schematic  depiction  of  the  opti¬ 
cal  transitions  used  for  spin  initialization,  control,  and 
readout.  These  optical  transitions  are  “indirect”  in  that 
they  couple  to  X+  states  with  an  electron  in  the  top  QD 
and  a  hole  in  the  bottom  QD.  These  indirect  transitions 
are  responsible  for  the  enhanced  wavelength  tunability 


we  achieve  with  this  qubit  design  and  control  strategy. 
In  Sect.  Ill  Alwe  analyze  the  tunability  and  optical  dipole 
strength  of  indirect  transitions.  In  Sect.  Ill  51  we  analyze 
the  spin  initialization  and  readout  protocols.  We  develop 
the  coherent  control  strategy  in  Sect.  Ill  Cl 


A.  Indirect  transitions 

Direct  transitions,  as  schematically  indicated  in 
Fig.  0}i,  involve  electrons  and  holes  in  the  same  QDs. 

These  direct  transitions  (e.g.  and  h&ve  a 

weak  dependence  of  their  energy  on  the  applied  electric 
field  due  to  the  quantum  confined  Stark  shift.  The  en¬ 
ergy  of  the  direct  transition  of  the  neutral  exciton  (A°, 
one  electron  and  one  hole)  in  a  QDM  with  QDs  sepa¬ 
rated  by  a  4  nrn  barrier  is  shown  by  the  black  symbols 
in  Fig.  [3lr.  Indirect  transitions,  also  depicted  in  Fig.  [l}\, 
involve  electrons  and  holes  in  separate  QDs.  Because  the 
relative  energy  levels  of  the  two  QDs  shift  in  response  to 

an  applied  electric  field,  indirect  transitions  (e.g. 

and  )  have  a  wavelength  that  depends  strongly  on 

the  applied  electric  field.  The  energy  of  PL  emitted  by 
an  indirect  X°  transition  is  shown  by  the  red  symbols 
in  Fig.  [3jc.  Fig.  [3jr  demonstrates  that  the  energy  of  an 
indirect  transition  can  be  tuned  by  at  least  18  meV  as 
the  applied  electric  field  is  varied  from  5  to  45  kV /cm. 

The  sensitivity  of  indirect  transition  energies  to  applied 
electric  field  can  be  enhanced  by  placing  a  thicker  barrier 
between  the  QDs.  However,  increasing  the  thickness  of 
the  barrier  also  decreases  the  optical  dipole  matrix  cle¬ 
ment  of  the  indirect  transition  by  reducing  the  overlap 
between  the  wavefunctions  of  the  electron  and  hole.  To 
assess  the  relative  strength  of  the  optical  dipole  matrix  el¬ 
ement  for  indirect  transitions,  in  Fig.  [HJo  we  plot  the  ratio 
of  PL  intensity  emitted  by  indirect  and  direct  transitions. 
The  intensity  of  PL  emission  is  not  a  direct  measure  of  the 
optical  dipole  strength  because  the  competing  dynamics 
of  radiative  recombination  and  carrier  relaxation  into  and 
between  discrete  states  influences  the  probability  that 
charges  occupy  the  different  optically  excited  states.  For 
example,  when  the  indirect  transition  is  at  higher  energy 
than  the  direct  transition,  we  would  expect  that  nonra- 
diative  relaxation  of  the  hole  (e.g.  -A  )  would 

favor  emission  from  the  direct  exciton  and  reduce  the 
intensity  of  the  indirect  transition.  In  Fig.  [3jo  we  ob¬ 
serve  significant  fluctuations  in  the  relative  intensity  of 
the  PL  emitted  by  indirect  and  direct  transitions.  The 
discontinuity  near  the  anticrossing  of  the  direct  and  in¬ 
direct  transitions  (AC1)  is  partially  due  to  the  fact  that 
the  states  are  fully  molecular  within  this  range  of  electric 
fields  and  the  distinction  between  direct  and  indirect  does 
not  apply.  The  fluctuation  in  the  intensity  ratio  around 
AC2  arises  because  the  intensity  of  the  direct  transition 
decreases  due  to  an  anticrossing  with  an  excited  indirect 
transition  (not  shown). 
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FIG.  3:  (Color  Online)  a)  Energy  of  direct  and  indirect  PL 
from  the  neutral  exciton  in  a  single  QDM  as  a  function  of 
applied  electric  field,  b)  Ratio  of  PL  intensity  (indirect  /  di¬ 
rect)  as  a  function  of  applied  electric  field.  AC1  indicates  the 
region  in  which  the  direct  and  indirect  transitions  anticross. 
AC2  indicates  a  region  in  which  the  energy  and  intensity  of 
the  direct  transition  is  influenced  by  an  anticrossing  with  an 
excited  indirect  transition  (not  shown). 


The  PL  intensity  ratio  shown  in  Fig.  [3jo  allows  us  to 
estimate  that  the  dipole  matrix  element  for  indirect  tran¬ 
sitions  in  this  QDM  has  an  average  value  approximately 
0.4  times  the  dipole  matrix  element  of  the  direct  transi¬ 
tion,  as  shown  by  the  dashed  horizontal  line  in  Fig.  (3Jd. 
We  note  that  there  can  be  significant  fluctuations  be¬ 
tween  QDMs  in  the  relative  intensity  of  indirect  transi¬ 
tions.  The  particular  QDM  studied  in  Fig.[3]has  a  some¬ 
what  large  indirect  transition  intensity.  This  increased 
intensity  for  the  indirect  transition  is  likely  correlated  to 
a  weak  tunnel  coupling  in  this  QDM  that  inhibits  hole 
relaxation  between  QDs.  This  conclusion  is  supported 
by  the  observation  of  a  relatively  small  anticrossing  en¬ 
ergy.  We  take  0.2  as  a  reasonable  approximation  of  the 
ratio  of  dipole  matrix  elements  for  indirect  and  direct 
optical  transitions  in  QDMs  with  a  4  nm  barrier.  This 
value  is  depicted  by  the  solid  red  horizontal  line  in  Fig.jSJr 
and  is  used  in  our  our  device  design  and  QDM  modeling. 
Further  investigation  of  the  relationship  between  QDM 
structure  and  the  dipole  matrix  element  for  indirect  tran¬ 
sitions  can  guide  engineering  of  QDMs  to  enhance  the 


coupling  of  indirect  transitions  to  optical  fields. 


B.  Initialization  and  readout 


Our  spin  initialization  protocol  (Fig.  |4}i)  is  based  on 
optical  shelving.  To  initialize  our  qubit  in  the  hole  spin- 
up  state,  we  illuminate  the  QDM  with  <j~  polarized  light 
in  resonance  with  the  optical  transition  that  excites  from 

the  (™)  state  to  one  molecular-like  branch  of  the  ad¬ 
mixture  ±  This  transition  is  indicated  by 

the  pump  laser  transition  in  Fig.  |4}i-  Because  the  opti¬ 
cally  excited  state  is  indirect,  the  energy  of  this 

optical  transition  tunes  very  strongly  with  applied  elec¬ 
tric  field,  enabling  the  wavelength  tunability  we  discuss 

below.  The  admixture  ±  (o’fKj.)  possible  be¬ 

cause  hole  spin  mixing  allows  the  spin-flip  tunneling  that 
couples  the  hole  spin  down  in  the  bottom  QD  with  the 
hole  spin  up  in  the  top  QD.  This  mixing,  and  the  conse¬ 
quent  formation  of  a  molecular  state,  is  indicated  by  the 
+/-  joining  the  two  atomic- like  basis  states.  The  relative 
weight  of  the  two  basis  states  does  depend  on  the  applied 
electric  field,  as  we  discuss  further  below. 

The  optically  excited  state  ±  may  radia- 

tively  decay  into  the  state  by  emission  of  a  di¬ 
rect  exciton.  is  a  metastable  state  in  the  range 

of  electric  fields  considered  here  because  the  hole  in  the 
ground  state  of  the  bottom  QD  is  at  higher  energy  than 
the  ground  state  of  the  top  QD.  The  hole  will  thus  relax 
back  to  the  top  QD  and  be  again  subject  to  the  spin  ini¬ 
tialization  laser.  As  a  result  of  the  hole  spin  mixing  there 
is  approximately  1%  probability  for  radiative  relaxation 

from  to  the  (o’°)  ground  state.  After  such 

a  radiative  decay  the  hole  spin  state  has  been  initialized 
to  the  spin  up  projection.  The  narrow  band  pump  laser 
is  not  resonant  with  any  transitions  that  couple  to  the 
hole-spin-up  state,  so  the  hole  spin  is  rapidly  shelved  in 
the  spin-up  projection. 

Spin  readout  is  achieved  by  measuring  resonance  fluo¬ 
rescence  as  depicted  in  Fig.  0fc.  The  readout  states  are 
immune  to  hole  spin  mixing  because  the  Pauli  exclusion 

principle  forbids  any  state  with  form  The  read¬ 


out  states  depicted  (e.g.  )  are  m  fact  molecular 

states  (e.g.  ±  (o’tii))’  but  this  spin-conserving 


mixing  does  not  alter  the  polarization  selection  rules  or 
enable  spin  flips  of  the  logical  basis  states  during  read¬ 
out.  The  direct  optical  transitions  (electron  and  hole  in 
same  QD)  of  the  readout  states  are  strongly  suppressed 
because  the  parallel  projections  of  electron  and  hole  spin 
are  dark  exciton  configurations  that  do  not  couple  to  op¬ 
tical  fields.  The  only  mechanism  for  degradation  of  the 
readout  states  is  a  spin-flip  of  the  electron,  which  can  be 
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a)  Spin  Initialization 


b)  Coherent  Control 


c)  Nondestructive  Readout 


FIG.  4:  (Color  Online)  Proposed  strategy  for  a)  spin  initialization,  b)  manipulation,  and  c)  readout  using  optical  transitions  of 
QDMs.  Top  (bottom)  rows  indicate  the  spin  projections  of  electrons  (holes)  in  the  top  (right  column)  and  bottom  (left  column) 
QDs. 


energetically  suppressed  by  choosing  to  readout  from  the 
lower  energy  (ct+)  transition.  The  spin  readout  transi¬ 
tions  are  thus  cycling  transitions  compatible  with  nonde¬ 
structive  readout:  many  cycles  of  optical  absorption  and 
emission  can  be  undertaken  to  increase  detection  proba¬ 
bility  without  risk  of  altering  the  spin  projection  of  the 
logical  basis  stated 

Both  the  spin  initialization  and  spin  readout  protocols 
utilize  narrowband  lasers  resonant  with  specific  optical 
transitions.  Because  the  optical  transitions  used  for  both 
initialization  and  readout  are  indirect,  the  energy  of  this 
optical  transition  can  be  varied  with  the  applied  electric 
field  to  tune  the  optical  transitions  of  an  arbitrary  QDM 
into  resonance  with  narrow  band  initialization  and  read¬ 
out  lasers  utilized  for  multiple  QDMs  within  the  same 
device.  In  Fig.  [5ji  we  plot  the  energy  of  optical  transi¬ 
tions  in  the  QDM  designed  for  this  device  strategy.  The 
optical  transitions  are  calculated  by  taking  the  difference 
between  the  calculated  X+  and  h+  state  energies  plotted 
in  Fig.  [2]  We  plot  results  in  Fig.  [5] only  for  electric  fields 
above  Fp,  where  the  logical  basis  states  remain  suitable 
for  all-optical  control. 

The  results  plotted  in  Fig.  [5}a  demonstrate  that  the  in¬ 
direct  optical  transitions  can  be  tuned  by  3  meV  (from 
1285  to  1288  meV)  with  applied  electric  fields  ranging 
between  approximately  12  and  20  kV/cm.  The  linear  de¬ 
pendence  of  wavelength  on  applied  electric  field  continues 
far  beyond  20  kV/cm,  but  we  restrict  Fig.  [5]to  this  range 
of  electric  fields  for  clarity  in  the  discussion  of  the  coher¬ 
ent  control  protocol.  The  maximum  wavelength  tunabil- 
ity  that  can  be  achieved  with  indirect  transitions  used 
for  spin  initialization  and  readout  will  likely  be  limited 
by  the  applied  electric  field  at  which  the  lowest  confined 
hole  state  of  the  bottom  QD  crosses  the  first  excited  hole 
state  of  the  top  QD,  at  which  point  additional  spin  in¬ 
teractions  and/or  relaxation  mechanisms  could  become 
important.  Spectroscopy  of  the  excited  states  of  holes  in 
QDMs  suggests  that  this  limit  would  permit  tuning  over 


approximately  10  meV4£ 

Tuning  over  10  meV  is  an  order  of  magnitude  improve¬ 
ment  over  the  typical  Stark  shift  tuning  (of  order  1  meV) 
achieved  for  single  InAs  QDsi^  and  comparable  to  the  gi¬ 
ant  Stark  shift  that  can  be  achieved  for  single  QDs  con¬ 
fined  between  AlGaAs  barriers— i  In  our  QDM  design  the 
AlGaAs  layer  that  blocks  injection  of  carriers  from  the 
n-type  GaAs  can  be  moved  arbitrarily  far  away  from  the 
QDs.  As  a  result,  our  QDM  design  is  likely  to  be  more 
compatible  with  the  fabrication  of  photonic  crystal  cav¬ 
ities  than  single  QDs  that  achieve  a  giant  Stark  shift  by 
placing  AlGaAs  layers  in  close  proximity  to  the  QDs.  As 
discussed  in  Sect.  Ill  Al  the  dipole  matrix  element  for  in¬ 
direct  transitions  is  weaker  than  the  comparable  dipole 
matrix  element  for  a  direct  transition.  We  discuss  the 
impact  of  these  dipole  matrix  elements  in  more  detail 
below. 


C.  Spin  rotations 

In  this  section  we  develop  our  approach  to  the  im¬ 
plementation  of  single  qubit  rotations.  We  quantify  the 
quality  of  our  gates  by  calculating  the  fidelities  of  our 
spin  rotations  in  the  presence  of  loss  mechanisms  and 
unwanted  dynamics. 

To  implement  arbitrary  spin  rotations  we  must  be  able 
to  perform  rotations  about  two  orthogonal  axes;  rota¬ 
tions  about  other  axes  can  be  implemented  by  combining 
those.  In  our  present  proposed  design,  rotations  about 
the  quantization  (z)  axis  will  be  carried  out  using  the 
cycling  (measurement)  transitions:  a  cyclic  pulse  on  the 
1 7 ~  measurement  transition  will  induce  a  phase  and  im¬ 
plement  a  rotation  about  the  z  axis  by  an  angle  deter¬ 
mined  by  the  detuning.— £  The  detuning  can  be  controlled 
independently  for  each  qubit  state  (defined  by  an  indi¬ 
vidual  QDM)  using  the  electric  field  that  tunes  the  indi¬ 
rect  transition  relative  to  a  fixed  optical  source  frequency. 
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FIG.  5:  (Color  Online)  a)  Energies  of  optical  transitions  in 
our  designed  QDM  as  a  function  of  applied  electric  field.  The 
two  top  most  lines  (red)  indicate  the  transitions  between  the 
two  qubit  states  and  the  target  A'+  state  (| t}).  The  second 
highest  pair  of  lines  (black)  indicate  the  energy  of  transitions 
involving  the  unwanted  X+  state  (|u)).  Other  colors  indicate 
the  energies  of  all  other  transitions  with  the  same  polarization 
(<r-)  that  have  been  included  in  our  calculations.  Inset  indi¬ 
cates  where  these  transitions  lie  relative  to  the  anticrossing 
of  the  X+  state,  b)  Dipole  strengths  of  each  optical  tran¬ 
sition  as  a  function  of  applied  held.  Color  and  dashing  as 
in  top  panel.  Dipole  strengths  are  measured  relative  to  the 
dipole  for  a  direct  transition,  which  is  given  a  strength  of  1. 
c)  Calculated  fidelity  of  a  n/2  rotation  using  the  two  target 
transitions  in  the  presence  of  all  parasitic  transitions  and  tak¬ 
ing  into  account  dipole  variations.  The  lower  axis  indicates 
the  QDM’s  detuning  from  an  ideal  laser  for  each  value  of  the 
applied  electric  held. 


Coherently  exciting  the  two  target  (Raman)  transitions 
along  the  quantization  axis  with  circularly  polarized  light 
allows  us  to  implement  rotations  about  an  orthogonal 
axis. 

To  implement  rotations  about  the  x  axis  we  pro¬ 
pose  the  use  of  the  A  system  formed  by  the  two  qubit 

states  ( (o’”)  and  (”’!j|) )  and  one  °f  the  spin-mixed  trion 
states,  as  depicted  in  Fig.  0Jd.  As  shown  in  Fig.  HJo, 


there  are  two  excited  states  (the  spin-mixed  trion  states 
a±  (j}’{l)  P±  (o’-flqi.) )  that  are  very  close  in  energy.  We 
develop  our  strategy  based  on  the  use  of  the  high-energy 
spin-mixed  trion  state  for  the  Raman  transition  and  la¬ 
bel  this  state  1 1),  for  target.  We  want  to  avoid  the  low- 
energy  spin- mixed  trion  state  and  label  this  state  |u),  for 
unwanted.  We  will  return  to  the  impact  of  this  unwanted 
state  below. 

The  two  atomic-like  basis  states  that  contribute  to  the 
spin-mixed  molecular  state  |f)  involve  direct  and  indirect 
transitions  with  significantly  different  optical  dipoles.  A 
key  advantage  of  our  scheme  is  that  the  molecular  state 
that  mediates  the  Raman  transition  is  dominated  by  the 
triplet  state  that  is  excited  by  the  indirect  transition 
( (U))-  Hole  spin  mixing  adds  a  small  fraction  of  the 

singlet  state  excited  by  the  direct  transition  ( (o'^)  )•  Be¬ 
cause  the  molecular  state  is  dominated  by  the  triplet  con¬ 
figuration  basis  state,  the  dipole  matrix  element  for  the 
indirect  transition  is  comparable  to  the  direct  transition. 
Both  transitions  have  the  same  polarization.  The  en¬ 
ergy  levels  and  dipole  strengths  are  computed  assuming 
a  hole  spin  mixing  strength  of  300  /mV,  which  we  believe 
to  be  a  value  achievable  with  present  growth  methods 
(see  Sect.  El  and  Appendix  iBl) . 

Spin  rotations  about  the  x  axis  are  based  on  coherent 
population  trapping  (CPT).  As  developed  in  Ref.  @  we 
can  create  an  effective  two-level  system  by  using  two  ap¬ 
propriate  phase-locked  pulses,  each  focused  on  one  of  the 
legs  of  the  A  system.  We  create  an  effective  two-level  sys¬ 
tem  when  the  parameters  of  the  two  pulses  are  related  in 
a  certain  way£  physically,  a  dark  state  is  created.  When 
the  two  transitions  have  the  same  Rabi  frequency,  the 

bright  state  is  an  equal  superposition  of  (o’jj.)  and  (o  il)  ■ 
Therefore  a  cyclic  evolution  will  induce  a  phase  to  this 
superposition  state  and  implement  a  rotation  about  the 
x  axis.  The  angle  of  rotation  can  be  controlled  by  the 
detuning  {e  in  Fig.  (4jp) ,  which  is  the  same  for  both  legs. 

There  are  always  sources  of  error  that  will  lower  the 
fidelity  of  the  spin  rotations.  The  finite  lifetime  of  the 
excited  state  will  introduce  non  unitary  evolution  and 
therefore  will  reduce  the  purity  of  the  qubit.  The  un¬ 
wanted  state  |m),  along  with  other  parasitic  transitions, 
will  effect  unitary  but  unwanted  dynamics  and  will  also 
lower  the  quality  of  the  gate.  Finally,  there  is  an  ad¬ 
ditional  error  coming  from  the  fact  that  the  pulse  that 
is  meant  to  excite  one  transition  of  the  A  system  will 
also  affect  the  other  one.  This  is  because  the  two  transi¬ 
tions  have  the  same  polarization.  To  quantify  the  effects 
of  these  error  mechanisms  and  compare  our  approach  to 
the  ‘traditional’  approach  in  which  the  control  lasers  are 
tuned  relative  to  the  fixed  frequency  of  QD  optical  tran¬ 
sitions,  we  calculate  the  fidelity  of  the  gate.  The  fidelity 
is  calculated  as  (Tr(ptp))^,  where  pt  is  the  density  ma¬ 
trix  corresponding  to  the  target  final  state  when  starting 
from  initial  state  dt,  p  is  the  actual  final  density  ma¬ 
trix,  and  the  average  is  taken  over  all  initial  spin  states. 


FIG.  6:  Fidelity  as  function  of  angle  of  rotation  for  ‘tradi¬ 
tional’  approach  (dashed  black  line)  in  which  the  laser  is  de¬ 
tuned  and  our  new  approach  (solid  red  line)  in  which  the 
QDM  is  tuned  relative  to  a  fixed  laser.  The  7r-rotation  is  im¬ 
plemented  by  a  resonant  set  of  pulses,  in  which  case  the  two 
approaches  coincide.  Inset:  temporal  profile  of  the  hyperbolic 
secant  pulses  used  in  our  simulations  (intensity  in  a.u.).  The 
bandwidth  of  the  pulse  is  0.02  meV;  the  optical  decay  time  is 
taken  to  be  Ins. 


The  fidelity  for  the  ‘traditional’  approach  is  plotted  by 
the  dashed  black  line  in  Fig.  [G]  as  a  function  of  the  an¬ 
gle  of  the  rr-rotation.  The  angle  of  rotation  is  controlled 
by  varying  the  detuning  of  the  laser  while  the  electric 
field  applied  to  the  QDM,  and  thus  the  optical  transition 
energies  of  the  QDM,  remain  fixed.  We  choose  a  fixed 
electric  field  of  15.3  kV/cm,  where  the  dipole  strengths 
of  the  two  branches  of  the  Raman  transition  are  equal 
(see  Fig.  [SJd).  The  inset  to  Fig.  [G]  shows  the  temporal 
pulse  envelope  used  in  these  calculations.  In  the  next 
section  we  assess  the  impact  of  changing  dipole  strengths 
and  electric  fields  on  the  fidelity  of  tunable  spin  rotations 
implemented  using  our  approach. 

D.  Tuning  spin  rotations  for  a  single  QDM 

As  discussed  above,  the  angle  of  coherent  spin  rota¬ 
tions  can  be  controlled  by  varying  the  detuning  of  the 
lasers  relative  to  the  optical  transitions  (e  in  Fig.  |4]d)  .  In 
traditional  device  designs  the  laser  would  be  tuned  rel¬ 
ative  to  the  fixed  frequency  of  optical  transitions  of  the 
QD  or  QDM  that  defined  the  qubit.  The  fidelity  of  gate 
operations  that  can  be  achieved  for  our  QDMs  using  the 
traditional  approach  is  plotted  by  the  dashed  black  line 
in  Fig.  [6]  Our  design  achieves  improved  scalability  by 
utilizing  laser  sources  with  fixed  frequency  and  tuning  in¬ 
dividual  qubits  into  resonance  with  the  available  sources. 
The  tunability  can  be  used  both  to  select  which  qubits 
are  affected  by  the  lasers  and  to  control  the  detuning  in 
order  to  control  the  angle  of  rotation.  Unfortunately,  the 
applied  electric  field  that  tunes  the  energies  also  changes 
the  dipoles  of  the  transitions  because  the  composition  of 


the  states  varies.  As  a  result,  the  dipoles  do  not  remain 
constant  when  we  vary  the  applied  electric  field  in  order 
to  control  the  detuning.  These  varying  dipoles  are  shown 
in  Fig.  [5Jd. 

To  show  that  the  variation  in  dipole  matrix  elements 
does  not  prohibit  execution  of  spin  rotations  by  arbitrary 
angles,  we  calculate  the  fidelity  for  spin  rotations  exe¬ 
cuted  by  tuning  QDMs  relative  to  a  fixed  frequency  laser 
source.  The  computed  fidelities  are  shown  by  the  solid 
red  line  in  Fig.  [6]  The  calculations  include  all  sources  of 
error  considered  for  the  traditional  approach  in  addition 
to  the  changing  dipole  strengths.  The  laser  pulse  profile 
and  frequency  remain  constant  and  the  angle  of  rotation 
is  varied  by  tuning  the  electric  field  to  control  the  detun¬ 
ing  (e).  The  plotted  value  indicates  the  best  fidelity  that 
can  be  achieved  for  a  given  angle  of  rotation  using  the 
design  from  Ref.  @  without  additional  pulse  shaping.  To 
compensate  for  the  unequal  dipoles  we  vary  the  relative 
intensity  of  the  two  lasers  that  address  the  two  legs  of 
the  A  transition.  Rapid  modulation  of  laser  power  for 
each  target  rotation  will  be  challenging,  but  should  be 
achievable  with  electrooptic  modulators.  The  optical  de¬ 
cay  time  is  taken  to  be  constant  at  1  ns  because  direct 
radiative  recombination  will  likely  remain  the  dominant 
lifetime  limit.  Fig.[6]shows  that  the  fidelity  of  gate  opera¬ 
tions  implemented  with  our  approach  is,  in  fact,  slightly 
better  than  that  achieved  by  the  traditional  approach. 
The  electric  field  must  be  tuned  by  only  0.05  kV/cm  to 
affect  the  change  from  a  7r/2  to  a  tt  rotation.  Fig.  [louses 
simple  pulse  shapes  and  therefore  provides  a  conservative 
estimate  of  the  fidelities  that  can  be  achieved  with  our 
approach. 


III.  SCALABILITY 

In  this  section  we  address  two  important  aspects  of 
scalability.  First,  we  examine  the  tunability  of  the  sys¬ 
tem,  in  the  sense  of  fixing  a  modest  number  of  optical 
sources  and  tuning  the  QDMs  into  resonance  with  those 
sources.  This  approach  allows  for  a  substantial  decrease 
in  the  equipment  overhead  and  device  complexity.  Sec¬ 
ond,  we  develop  an  approach  for  coupling  arbitrary  pairs 
of  qubits  without  increasing  our  laser  overhead. 


A.  Tunability  in  a  multi-QDM  system 

In  the  ideal  case,  all  QD-based  qubits  would  be  iden¬ 
tical  and  could  be  addressed  by  a  single  group  of  optical 
sources.  In  reality,  the  inhomogeneous  distribution  of  en¬ 
ergy  levels  in  QDs  makes  this  impossible.  A  mechanism 
for  tuning  the  optical  initialization,  control,  and  readout 
transitions  into  resonance  with  a  small  number  of  fixed- 
frequency  optical  sources  is  therefore  desirable  for  the  de¬ 
velopment  of  a  realistic,  scalable,  multi-qubit  device.  We 
have  analyzed  the  wavelength  tunability  of  the  spin  ini¬ 
tialization  and  readout  transitions  above,  and  addressed 
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rotation  angle  (n) 


FIG.  7:  (Color  online)  Fidelity  of  spin  rotation  as  a  function 
of  the  angle  of  rotation  for  a  multi-QDM  system  with  varia¬ 
tions  in  their  optical  resonances.  Each  curve  corresponds  to  a 
different  value  of  the  E-field;  from  top  to  bottom  at  (f>  =  7t/2: 
13.8,  14.7,  15.3,  15.6,  16,  16.4,  16.8,  17.2,  13.4,  17.7,  18.1,  19.3 
kV/cm.  The  pulse  envelope  is  similar  to  that  used  in  Fig.  [6] 
but  with  half  the  bandwidth,  so  about  double  in  temporal 
length. 


spin  rotation  tunability  within  a  single  qubit.  We  now 
address  the  wavelength  tunability  of  our  spin  rotations 
in  the  case  of  a  system  of  inhomogeneous  QDMs. 

The  inhomogeneous  distribution  of  QD  energy  lev¬ 
els  means  that  an  arbitrary  QDM  qubit  will  not  have 
equal  dipoles  for  the  two  legs  of  the  A  transitions  when 
that  QDM  qubit  is  tuned  into  resonance  with  the  laser 
sources.  As  a  result,  Fig.  [6]  is  not  general  enough  to  de¬ 
scribe  a  multi-QDM  system.  To  analyze  the  multi-QDM 
system  we  compute  the  fidelity  versus  rotation  angle  for 
various  values  of  the  E- field,  which  simulates  the  need 
to  tune  an  individual  QDM  farther  away  from  the  ideal 
electric  field  (at  which  the  dipole  strengths  of  the  two 
transitions  are  equal)  in  order  to  tune  that  QDM  into 
resonance  with  the  available  optical  source.  The  results 
are  plotted  in  Fig.  0  The  different  lines  correspond  to 
different  applied  electric  fields  with  the  consequent  differ¬ 
ent  relative  dipole  strengths  for  the  A  transition.  In  this 
calculation  we  do  not  correlate  the  detuning  of  the  tar¬ 
get  transitions  to  the  dipoles.  Fig.  [7]  therefore  describes 
the  range  of  fidelities  that  could  be  achieved  for  a  QDM 
ensemble  with  a  certain  energy  distribution. 

To  further  quantify  the  range  of  tuning  that  could  be 
achieved  for  an  arbitrary  QDM  qubit,  we  plot  the  fidelity 
for  a  7t/2  rotation  as  a  function  of  electric  field  applied  to 
tune  the  QDM  into  resonance  with  a  fixed  laser  source. 
The  results,  presented  in  Fig.  [5jy  show  that  our  scheme 
can  achieve  a  fidelity  of  at  least  0.8  when  tuning  a  partic¬ 
ular  QDM  to  use  a  fixed-frequency  laser  anywhere  within 
a  2.5  meV  window.  A  fidelity  of  0.8  is  certainly  well  be¬ 
low  the  formidable  values  (~99.99%)  needed  for  practical 
quantum  computation.  The  present  fidelity  values  are 
intended  to  illustrate  that  is  feasible  to  tunc  individual 


QDMs  to  utilize  fixed-frequency  laser  sources  and  should 
be  viewed  as  a  lower  bound  of  what  can  be  achieved.  The 
present  calculations  have  used  very  simple  pulse  shapes 
and  there  are  ways  to  increase  the  fidelity  of  the  oper¬ 
ations  by  using  numerical  pulse  optimization  techniques 
to  design  pulse  shapes  tailored  to  this  system.  In  partic¬ 
ular,  by  use  of  Optimal  Control  Theory  we  can  impose 
the  condition  that  the  rotation  is  robust  against  substan¬ 
tial  variations  in  the  ratio  of  the  dipole  couplings.  This 
could  be  done  by  running  a  Krotov  algorithm  in  parallel 
for  a  number  of  systems  with  varying  dipoles.  Similar 
approaches  can  improve  fidelities  against  other  types  of 
errors,  such  as  unintended  couplings  to  other  states.  We 
expect  that  with  such  techniques  the  fidelities  should  ap¬ 
proach  the  near-perfect  values  needed  for  realistic  quan¬ 
tum  information  processing. 


B.  Two-qubit  gates 

The  wavelength  tunability  enabled  by  our  QDM-based 
qubit  design  provides  a  unique  opportunity  to  couple  ar¬ 
bitrary  pairs  of  qubits  through  a  photonic  crystal  cavity 
mode.  As  depicted  in  F ig . |TJ:  and  d,  the  qubits  in  our  pro¬ 
posed  device  architecture  are  defined  by  QDMs  that  are 
deterministically  spatially-coupled  to  a  single  photonic 
crystal  membrane  cavity  models  The  wavelength  of  the 
photonic  crystal  mode  can  be  fine  tuned^  but  the  QD 
inhomogeneity  makes  it  extremely  challenging  to  choose 
a  cavity  resonance  that  is  ideal  for  all  qubits  if  they  are 
defined  by  as-grown  QDs  or  QDMs.  Our  QDM  qubit  de¬ 
sign  allows  each  QDM  to  be  individually  controlled  by 
local  electric  field  gates  to  provide  in  situ  tuning  of  the 
coupling  between  individual  QDMs  and  the  cavity  mode. 
We  propose  to  use  p-  and  n-type  GaAs  as  the  electrical 
contacts  in  order  to  avoid  the  inclusion  of  metals,  which 
would  heavily  degrade  the  cavity  quality  factor:-44  Dur¬ 
ing  device  calibration,  the  local  electric  field  necessary  to 
tune  each  QDM  into  resonance  with  the  available  optical 
sources  and/or  cavity  mode  would  be  determined.  Dur¬ 
ing  the  execution  of  logic  operations,  each  QDM  can  be 
tuned  appropriately  with  respect  to  these  optical  sources 
and  the  cavity^ 

Using  the  approach  of  Ref.  (HI,  pairwise  entanglement 
can  be  created  via  cavity-mediated  entangling  gates  by 
selectively  tuning  two  QDMs  to  the  appropriate  detuning 
relative  to  the  cavity  mode.  The  remaining  qubits  are  iso¬ 
lated  by  significantly  detuning  them  relative  to  the  cavity 
via  use  of  their  local  electric  fields.  The  cavity-mediated 
entanglement,  described  in  more  detail  below,  occurs 
only  when  a  laser  pulse  is  present,  suppressing  the  pos¬ 
sibility  of  unintended  couplings  as  the  QDMs  are  tuned 
in  preparation  for  each  logical  gate  operation.  The  opti¬ 
cal  transitions  relevant  to  coherent  single-spin  rotations, 
readout,  and  two-qubit  gates  are  all  among  the  highest- 
energy  transitions  shown  in  Fig.  [SJi.  Consequently,  the 
QDMs  can  be  red-detuned  relative  to  the  cavity  when 
they  are  to  be  isolated,  ensuring  that  other  optical  tran- 
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sitions  of  the  QDM  never  become  more  strongly  coupled 
to  the  cavity  mode. 

Our  proposed  cavity  mediated  two-qubit  gates  are 
based  on  the  readout  transitions  -H-  ± 

(o’-fHI-)  ■  Bringing  these  transitions  near  the  two-photon 

resonance  with  the  cavity^  we  can  diagonalize  the 
Hamiltonian  consisting  of  one  qubit  state  and  the  direct 
trion  excited  by  a  fixed  circular  polarization  for  each  of 
two  QDMs  and  by  three  photonic  number  states:  the 
zero,  one  and  two  photon  states.  The  eigenstates  of 
this  Hamiltonian  will  be  entangled  states  of  all  three 
systems  (QDM1,  QDM2,  photon),  and  will  naturally 
separate  into  subspaces,  depending  on  the  number  of 
excitations—  These  eigenstates  have  energies  that  are 
shifted  from  the  sum  of  the  non-interacting  system’s  net 
energy,  and  thus  a  specific  entangled  state  can  be  selec¬ 
tively  excited  optically.  The  requirement  on  the  cavity 
is  that  its  Q  is  high  enough  that  its  line  broadening  is 
smaller  than  the  Zeeman  energy  of  the  qubits  and  that 
the  QD-cavity  coupling  g  is  larger  than  the  decay  rates 
involved.  These  requirements  ensure  that  the  other  qubit 
states  will  not  couple  to  the  cavity  and  each  other  and 
that  two-qubit  unitary  operations  are  in  principle  possi¬ 
ble.  Under  these  conditions,  we  can  use  a  single  cyclic 
pulse  to  induce  a  tt  phase  on  the  state  |  tl'itl'2),  by  which 
we  mean  the  state  of  two  qubits,  each  in  their  own  QDM, 
both  with  hole  spin  up.  Inducing  a  tt  phase  on  the  state 
|  fhlta)  effects  an  entangling  CZ  gate— 


IV.  EXPERIMENTAL  EVIDENCE  OF 
HOLE-SPIN  MIXING  IN  OPTICALLY  EXCITED 
STATES  OF  QDMS 

As  a  first  step  toward  implementing  this  device  ar¬ 
chitecture,  we  experimentally  measure  the  magnitude  of 
hole  spin  mixing  in  the  positive  trion  (X+)  state.  The 
QDM  studied  experimentally  has  a  4  nm  barrier  sepa¬ 
rating  the  two  QDs,  identical  to  the  barrier  in  the  QDM 
proposed  here.  The  QDM  we  measured  is  grown  on  a  n- 
type  substrate  and  is  optically  (not  electrically)  charged 
with  a  single  excess  hole— ^  The  built-in  electric  field 
and  QD  asymmetry  of  this  sample  are  also  inverted  rela¬ 
tive  to  the  sample  described  above.  None  of  these  changes 
impact  the  spin  mixing  interactions  that  are  the  focus  of 
this  experiment.  The  sample  is  prepared  with  an  ohmic 
back  contact  and  Schottky  top  contacts  enabling  us  to 
study  individual  QDMs  using  photoluminescence  (PL) 
spectroscopy  as  a  function  of  both  applied  electric  and 
magnetic  fields.  A  detailed  description  of  experimental 
procedures  can  be  found  in  previous  publications—  Here 
we  present  the  evidence  of  hole-spin  mixing  in  the  opti¬ 
cally  excited  trion  states  (A+),  which  manifests  as  new 
anticrossings  in  the  photoluminescence  (PL)  spectra.  We 
compare  computational  simulations  of  the  observed  PL 
spectra  to  extract  a  numerical  value  of  the  hole  spin  mix¬ 
ing  strength. 


In  Fig.  [HJo  we  show  the  experimentally  measured  PL 
spectra  of  a  single  QDM  in  the  presence  of  a  6  T  mag¬ 
netic  field.  The  figure  plots  the  photoluminescence  in¬ 
tensity  (grayscale)  as  a  function  of  both  applied  electric 
field  (x  axis)  and  energy  (y  axis).  The  initial  and  fi¬ 
nal  states  of  optical  recombination  are  the  X+  and  h+ 
states,  analogous  to  those  displayed  in  Fig.  [5Jr  and  b. 
As  seen  in  Fig.  [2j  Coulomb  interactions  cause  the  anti¬ 
crossings  of  the  X+  and  h+  states  to  occur  at  different 
values  of  the  applied  electric  field.  As  a  result,  the  PL 
energy  as  a  function  of  electric  field  shows  the  charac¬ 
teristic  “x”  shaped  pattern—  Note  that  the  “x”  appears 
doubled  because  the  magnetic  field  introduces  a  Zeeman 
splitting  between  states  that  have  a  total  flip  of  all  spins. 
The  circled  regions  indicate  the  locations  of  anticrossings 
and  fine  structure  that  arise  due  to  hole  spin  mixing. 


FIG.  8:  (Color  Online)  Experimental  and  calculated  spectra 
from  the  positive  trion  in  a  QDM  with  a  4  nm  barrier  in  a 
6  T  longitudinal  magnetic  field,  a)  Experimental  spectra,  b) 
Calculation  with  hm  =  0,  c)  Calculation  with  hm  =  0.2  meV. 
The  circles  highlight  four  regions  where  the  features  observed 
experimentally  are  only  reproduced  by  the  calculation  includ¬ 
ing  nonzero  hm. 

In  Figs.  [HJo  and  c  we  show  calculated  PL  spectra  com¬ 
puted  with  the  matrix  Hamiltonian  method  described  in 
Appendix  [B]  and  in  previous  publications—  In  Fig.  |5Jd 
the  hole  spin  mixing  term  (hm)  is  set  to  zero  and  in 
Fig.  [HJc  the  hole  spin  mixing  term  is  set  to  0.2  meV.  It 
is  clear  that  the  inclusion  of  the  hole  spin  mixing  term 
generates  the  additional  anticrossings  and  fine  structure 
apparent  in  the  circled  regions.  The  value  of  hm  ~  0.2 
meV  is  obtained  by  varying  the  magnitude  of  hm  and 
evaluating  the  best  fit  to  the  experimental  data.  This 
result  experimentally  verifies  the  existence  of  hole  spin 
mixing  in  the  trion  state  and  validates  the  feasibility  of 
the  qubit  design  and  device  architecture  proposed  here. 

V.  CONTROLLABLY  GENERATING  HOLE 
SPIN  MIXING  WITH  THE  DESIRED 
MAGNITUDE 

For  the  qubit  and  device  designs  we  propose,  strong 
spin  mixing  is  desirable  because  it  separates  the  target 
(| t))  and  unwanted  (|u))  states.  The  energy  level  cal¬ 
culations  presented  in  Fig.  [2]  use  a  value  of  hm  =  0.3 
meV,  and  the  energy  levels  and  dipole  matrix  elements 
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computed  with  these  parameters  are  the  basis  for  the  cal¬ 
culated  fidelity  presented  in  Fig.  [5]:.  The  fidelity  results 
would  improve  for  larger  values  of  hm. 

As  discussed  in  Ref.  [H,  holes  are  described  as  Lut- 
tinger  spinors  with  coupled  heavy  hole  (HH)  and  light 
hole  (LH)  components.  The  hole  pseudospin  projections 
in  self-assembled  QDs,  ff  and  4J-,  can  be  identified  with 
the  total  angular  momentum  of  the  Luttinger  spinors, 
Fz  =  3/2  and  Fz  =  —3/2.  When  the  axial  symmetry 
of  the  nanostructure  is  broken,  Fz  is  no  longer  a  good 
quantum  number  and  hole  spin  mixing  (i.e.  coupling 
between  Fz  =  3/2  and  Fz  =  —3/2  states)  arises.  The 
admixing  of  the  two  pseudospins  is  mediated  by  the  LH 
components  of  the  spinor.35  In  QDMs,  the  constituent 
QDs  are  usually  symmetric  enough  to  neglect  this  spin 
mixing  mechanism.  If  the  QDM  is  misaligned,  however, 
the  application  of  an  electric  field  forcing  the  hole  to  de¬ 
localize  over  the  two  QDs  leads  to  a  severely  asymmetric 
orbital  and  spin  mixing  becomes  important.  In  Fig.  [T|i 
we  show  a  cross-sectional  STM  of  stacked  InAs  QDs  that 
shows  the  possibility  of  symmetry  breaking  due  to  lat¬ 
eral  offset  along  the  stacking  axis.  Clearly,  the  larger  the 
lateral  offset  between  the  QDs  of  the  QDM,  the  stronger 
the  asymmetry  and  the  stronger  the  spin  mixing.  Using 
k-p  calculations  we  previously  showed  that  a  QDM  could 
give  rise  to  hm  —  0.1  meV  for  center-to-center  lateral 
offsets  of  ~  5  nmi£  A  5  nm  offset  is  unusually  large  but 
within  the  observed  ranged 

As  described  in  Sect.  El  we  have  now  experimentally 
measured  values  of  hm  up  to  ~  0.2  meV.  In  order  to 
further  increase  hm  one  could  use  larger  offsets,  but  this 
structural  parameter  is  unfortunately  difficult  to  control. 
Moreover,  increasing  the  offset  has  the  negative  side  ef¬ 
fect  of  reducing  the  tunneling  rate.35  In  this  section  we 
propose  a  more  efficient  way  to  enhance  the  spin  mixing 
by  using  QDs  with  larger  aspect  ratio  (i.e.  larger  ratio 
of  height  to  radius).  Larger  aspect  ratios  qualitatively 
increase  the  hole  spin  mixing  because  the  anisotropic 
masses  of  HHs  and  LHs  cause  the  LH  character  to  in¬ 
crease  with  increasing  aspect  ratio  in  QDs  grown  along 
the  [001]  (z)  direction.  As  a  result  of  the  increased  LH 
character  the  pseudospin  coupling  increases.  We  analyze 
the  strength  of  hole  spin  mixing  as  a  function  of  QDM 
structure  using  k-p  theory  and  demonstrate  that  this  ap¬ 
proach  provides  a  feasible  path  to  fabrication  of  QDMs 
with  hm  >0.3  meV. 

The  spin  mixing  parameter  hm  is  defined  as  the  ma¬ 
trix  element  coupling  states  with  holes  localized  in  oppo¬ 
site  QDs  and  with  opposite  spin.  We  analyze  the  depen¬ 
dence  of  this  spin  mixing  parameter  on  QDM  structure 
and  symmetry  with  a  simplified  Hamiltonian.  In  a  basis 

formed  by  |a)  =  (^’o)  and  | b)  =  (o’jj.)  ,  the  spin  mixing 
Hamiltonian  reads: 


/  Ea  +  edF  hm  \ 

y  hm  Et  J  ’ 


where  Ea  and  Et  are  the  energies  of  \a)  and  | b)  at  zero 
electric  field.  At  the  resonant  electric  field  (Ea  +  edF)  = 
Et  =  E0  and  the  two  hybridized  states  have  energies 
E±  =  Eq  ±  hm.  Thus,  hm  =  A/2,  where  A  =  E+  —  E- 
is  the  magnitude  of  the  spin  anticrossing  gap  between 

(r.)  (£»)• 

To  determine  A,  we  run  numerical  simulations  using 
the  same  theoretical  model  and  material  parameters  as  in 
Ref.  HH.  A  typical  hole  energy  spectrum  as  a  function  of 
the  external  electric  field  is  presented  in  Fig.  [Dfa).  The 
spectrum  is  calculated  with  a  magnetic  field  B  =  6  T 
to  separate  the  spin-conserving  anticrossings  ( F  ~  —2 
kV /cm)  from  the  spin  anticrossings  we  are  interested  in 
(F  ~  0  and  F  ~  —4  kV/cm,  see  arrows).  We  note  that  in 
our  formulation  of  the  k-p  Hamiltonian  there  are  no  off- 
diagonal  magnetic  terms,  so  the  strength  of  the  HH-LH 
coupling  terms  does  not  depend  on  54^ 

Fig.  (Da)  reveals  that  there  are  two  spin  anticrossings, 
Ao  and  Ai.  In  general  the  magnitude  of  the  anticrossings 
is  asymmetric,  with  Ao  >  Ai.  This  is  because  the  spin 
mixing  is  mediated  by  LHs  delocalized  over  the  entire 
QDM.  HHs  come  into  resonance  at  F  ~  —2  kV/cm,  but 
LH  states  come  into  resonance  at  electric  fields  closer  to 
0  kV/cm  because  of  their  lighter  vertical  mass.  Conse¬ 
quently,  the  Ao  anticrossing  is  closer  to  the  resonant  field 
of  LHs,  where  LH  delocalization  and  hence  spin  mixing 
are  maximized. 

In  Fig.  EKb)  we  plot  the  magnitude  of  Ao  for  different 
QDMs  as  a  function  of  the  lateral  offset.  The  QDMs  are 
made  of  lens-shaped  QDs  with  ff  =  2  nm  height,  radius 
R  and  interdot  barrier  D.  Because  we  are  interested  in 
strong  spin  mixing,  we  focus  on  large  (but  experimentally 
accessible)  offsets.  For  QDs  with  R  =  15  nm,  Ao  reaches 
maximum  values  of  0.2  meV  ( hm  =  0.1  meV),  which 
is  the  case  reported  in  Ref.  |35|.  For  QDs  with  R  =  13 
nm  Ao  is  systematically  larger,  and  hm  =  0.15  meV  is 
attained  for  a  5  nm  offset.  Last,  for  QDs  with  R  =  10 
nm  the  magnitude  of  the  spin  mixing  is  even  stronger 
and  Ao  =  0.6  meV  (hm  =  0.3  meV)  is  already  within 
reach. 

Similar  trends  to  those  reported  in  Fig.  [9][b)  are  ob¬ 
tained  for  Ai,  albeit  with  smaller  magnitude  (Ai  < 
A0/2).  One  could  think  of  increasing  the  QD  aspect 
ratio  by  growing  taller  QDs  instead  of  reducing  R.  How¬ 
ever,  this  also  reduces  the  tunneling  rate  of  the  QDM,  so 
reducing  the  lateral  size  is  preferable.  The  lateral  size  of 
InAs  QDs  can  be  controlled  experimentally  by  growing  on 
templated  surfaces^  It  may  also  be  possible  to  use  new 
methods  for  the  growth  of  QDMs  with  precisely  tailored 
three-dimensional  configurations  to  engineer  physical  in¬ 
teractions  that  maximize  wavelength  tunability  and  gate 
operation  fidelity^ 
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F  (kV/cm) 


lateral  offset  (Ang) 


FIG.  9:  (Color  Online)  (a)  Single  hole  energy  as  a  function 
of  the  electric  field  in  a  QDM  subject  to  B  =  6  T.  The  QDs 
have  H  =  2  nm,  R  —  15  nm,  D  =  1.8  nm  and  are  offset  by  4.2 
nm.  (b)  Spin  anticrossing  gap  as  a  function  of  the  QDM  offset 
for  QDs  with  different  aspect  ratio.  Triangles  correspond  to 
(. R,D )  =  (15,1.8)  nm,  circles  to  (R,  D)  =  (13,1.7)  nm  and 
squares  to  (R,  D)  =  (10, 1.6)  nm. 


VI.  SUMMARY 

We  have  presented  a  scalable  qubit  design  and  device 
architecture  based  on  the  spin  states  of  single  holes  con¬ 
fined  in  a  QDM.  The  QDM  qubit  enables  a  new  strat¬ 
egy  for  all-optical  coherent  control  via  indirect  optical 
transitions.  The  wavelength  tunability  of  indirect  transi¬ 
tions  allows  us  to  develop  a  device  architecture  in  which 
locally-applied  electric  fields  are  used  to  tune  individual 
QDMs  into  resonance  with  a  single  photonic  cavity  mode 
or  a  small  number  of  optical  sources  at  fixed  wavelengths. 
Thus,  our  design  offers  a  tremendous  reduction  in  the 
required  overhead  as  compared  to  schemes  based  on  sin¬ 
gle  QDs.  We  take  advantage  of  hole  spin  mixing  in  or¬ 
der  to  execute  spin  initialization  and  coherent  spin  rota¬ 
tions  with  indirect  optical  transitions  and  without  trans¬ 
verse  magnetic  fields.  The  absence  of  transverse  magnetic 
fields  enables  nondestructive  readout.  We  show  how  cy¬ 
cling  transitions  can  be  used  for  nondestructive  readout 
by  optically  exciting  trion  states  that  are  immune  to  hole 
spin  mixing.  We  demonstrate  that  ultrafast  coherent  spin 
rotations  can  be  achieved  by  tuning  the  indirect  transi¬ 


tions  of  individual  QDMs  into  resonance  with  a  single 
optical  pulse  or  a  pair  of  pulses.  We  show  that  these 
spin  rotations  have  a  minimum  fidelity  of  0.8  over  a  min¬ 
imum  tuning  range  of  2.5  meV,  despite  fluctuating  dipole 
strengths,  and  describe  a  path  to  the  development  of  op¬ 
tical  control  protocols  that  can  enhance  the  fidelity  and 
the  wavelength  tunability  of  this  approach.  We  demon¬ 
strate  experimental  and  theoretical  progress  toward  the 
development  of  QDM  materials  to  implement  the  designs 
presented  and  validated  here. 
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Appendix  A:  QDM  Design  parameters 

The  substrate  doping,  choice  of  QD  asymmetry  and  ap¬ 
plied  electric  field  can  be  used  to  control  the  total  charge 
occupancy  and  spatial  location  of  charges  in  a  QDM.~ 
The  qubit  we  present  here  is  designed  to  induce  controlled 
coupling  of  hole  energy  levels  in  a  range  of  electric  fields 
in  which  the  QDM  remains  charged  with  only  one  hole. 
The  thickness  of  the  undoped  GaAs  separating  the  QDs 
from  the  p-doped  region  (25  nm)  is  set  by  the  Coulomb 
blockade  strength  and  allows  for  deterministic  charging 
of  the  QDM  with  a  single  hole.  The  QD  asymmetry  (bot¬ 
tom  QD  truncated  to  2.8  nm,  top  QD  truncated  to  3.2 
nm)  is  chosen  to  ensure  that  hole  levels  come  into  reso¬ 
nance  for  the  sign  and  range  of  electric  fields  applied  in 
the  p-i-n  diode  structure.  This  asymmetry  and  sign  of 
the  applied  electric  field  also  ensure  that  the  electron  is 
located  in  the  top  QD. 

The  thickness  of  the  GaAs  separating  the  top  QD  from 
the  AlGaAs  blocking  layer  is  less  critical;  it  is  only  neces¬ 
sary  to  ensure  that  that  the  AlGaAs  layer  is  sufficiently 
thick  to  suppress  tunneling  of  electrons  from  the  n-doped 
GaAs  top  contact  and  sufficiently  far  away  from  the  QD 
to  avoid  altering  growth  dynamics  or  hole  confinement. 
Because  the  thickness  of  the  GaAs  and  AlGaAs  layers 
above  the  top  QD  are  flexible,  the  net  thickness  of  the 
heterostructure  can  be  chosen  to  provide  the  optical  con¬ 
finement  necessary  to  fabricate  a  high-Q  photonic  crys¬ 
tal  that  places  the  QDMs  at  electric  field  maxima.  The 
ARGan-j)  As  (x  >  0.7)  sacrificial  layer  is  included  below 
the  p-doped  GaAs  region  in  order  to  facilitate  an  under¬ 
cut  etch  that  releases  the  photonic  crystal  membrane. 
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Appendix  B:  Computational  Model 


The  energy  levels,  optical  transitions,  and  op¬ 
tical  dipoles  presented  above  are  calculated  using 
matrix  Hamiltonian  methods  presented  in  previous 
publications ^c21i32i35~—  The  relevant  basis  states  for  the 
optical  ground  state  (logical  basis  states)  are: 


/  o  o  \  /o  o\  /  o  o  \ 
[it  o) 


(Bl) 


In  this  basis,  the  matrix  Hamiltonian  that  describes 
the  hole  states  is  given  by: 


J 


/  ^g.BBgh  0 

0  -\gLBBgh 

- th  +  \gBBgb  hm 

y  hlTl  th  ^  g*B  B gb 


-th  +  \^BBgb 
hm 

dF  +  \gBBgh 

0 


—hm 

- th  ~  \ gBBgb 

0 


dF  -  \ gBBgh  j 


(B2) 


The  definitions  and  numerical  values  for  each  symbol  In  this  basis,  the  matrix  Hamiltonian  that  describes 
are  provided  in  Table  Q]  The  basis  states  for  the  positive  the  positive  trion  state  is  given  by: 
trion  (A+)  matrix  Hamiltonian  are: 


(  0  t  Wo  t\  (  o  (  0  t  W  0  f\ 
\o  u)  { n-  n )  s  v^^/r  V  ^  ^  J  \ ft  ft/ 

(B3) 


J 


Ex+ 


/  T  -  \gBBge 
—tx+ 

0 

hm 
\  hm 


- tX+ 


dF  -  \g,BBge 
J 
0 
0 


0 

J 

dF  -  \g,BBge 
0 
0 


Note  that  the  full  matrix  is  block  diagonal  with  a  sub¬ 
matrix  for  the  electron  spin  —1/2  that  exactly  parallels 
the  +1/2  case  presented  here.  The  parameter  definitions 
and  numerical  values  are  given  in  Table  [[] 

The  energies  of  hole  and  trion  states  are  computed 
by  calculating  the  eigenvalues  of  the  matrix  Hamiltonian 
as  a  function  of  electric  field,  F.  The  hole  spin  purity 
(Fig.  [2jc)  is  calculated  by  plotting  the  relative  contri¬ 
butions  of  each  basis  state  to  the  eigenvector  that  de¬ 
scribes  the  admixture  of  all  states  for  a  given  value  of  F. 
The  dipole  matrix  elements  are  computed  by  the  prod¬ 
uct  'bh  ■  S  ■  where  ’1+  and  Tx+  are  the  eigenvectors 
describing  the  hole  and  trion  states  at  a  given  value  of 


hm 

0 

0 


dF  +  J  -  ^gBS  (ge  +  2 gh) 


hm 
0 
0 
0 

dF  -  J  -  (g, 


\ 


2gh)J 

m. 


F.  S'  is  a  matrix  describing  the  selection  rules  for  all 
optical  transitions.  S  assigns  amplitude  1  to  all  bright 
direct  transitions  (i.e.  electron  and  hole  in  the  same  QD 
with  opposite  spin  projections),  amplitude  0  to  all  dark 
transitions  and  amplitude  ^  to  all  transitions  involving 
singlets  or  triplets  that  are  superpositions  of  a  bright  and 
a  dark  direct  transition.  Bright  indirect  transitions  in  S 
are  given  amplitude  0.2,  consistent  with  the  measured  ra¬ 
tio  of  optical  intensities  for  direct  and  indirect  transitions 
in  samples  where  the  QDs  are  separated  by  a  4  nm  bar¬ 
rier.  Indirect  transitions  that  involve  singlets  or  triplets 
that  are  superpositions  of  a  bright  and  a  dark  state  are 
given  amplitude  0.1. 
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Ex+ 

1280  meV 

Energy  of  the  trion 

-th 

-0.3  meV 

bare  hole  tunneling  strength 

—tx+ 

-0.45  meV 

hole  tunneling  strength  in  trion 
state 

d 

4  nm 

QD  separation 

F 

varies 

applied  electric  field 

hm 

0.3  meV 

hole  spin  mixing  strength 

Hb 

57.9  /reV/T 

Bohr  magneton 

B 

1  T 

Magnetic  Field 

9h 

-1.555 

hole  g  factor 

9b 

0.4 

Barrier  contribution  to  hole  g  factor 

9e 

-0.64 

electron  g  factor 

r 

3.2145  meV 

Energy  shift  due  to  Coulomb  inter¬ 
actions  when  all  three  charges  are 
in  the  same  QD 

J 

0.116  meV 

Electron-hole  exchange  energy 

TABLE  I:  Definitions  and  numerical  values  for  computational 
parameters. 
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